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Abstract 

5h . 

Various nonleptonic decay channels mediated by the quark-level subprocess b — > sss show hints of 
deviation from the Standard Model expectations. We analyse the double-vector decay B — > (j>K* with 
different generic new physics structures and find the constraints on the parameter spaces of new physics. 
The allowed parameter spaces are compatible with, but further narrowed down from, those obtained from a 
similar analysis using pseudoscalar modes. We also discuss further predictions for this channel as well as for 
B s —¥ 4>4>, and show how even a partial measurement of the observables may discriminate between different 
models of new physics. 
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It has been felt for a long time that the data from the nonleptonic decays of the B meson, mediated by the 
Qh. transition b — > sss, is not what one expects from the Standard Model (SM) with the Cabibbo-Kobayashi- 
Maskawa (CKM) paradigm for CP violation. On one hand this led to a vigorous exercise for the understanding 
of low-energy QCD dynamics, including different models for calculating the long-distance part of the amplitude, 
and on the other hand this made such nonleptonic modes an ideal testing ground for indirect signals of New 
Physics (NP) [TJ. 
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Let us make it clear here and now that no data, taken alone, is a clear indication for NP. Some of the discrepancies 
(like the value of sin(2/3) extracted from B — > 4>Ks) are compatible with the SM predictions at less than 2a, 
and some (like the abnormally large branching ratio of B — > rj ' K) may be explained by SM dynamics yet to be 
fully understood. However, all data, considered together, have a significant pull away from the SM predictions, 
and one may hope that this will emerge as an indirect signal of NP. This is more so since the b — ► sss transition 
involves quarks of the second and the third generations where the SM is less well-tested. 

In an earlier publication we have shown that from a model-independent analysis of B — > P1-P2 (two pseu- 
doscalars) and B — > PV (one pseudoscalar and one vector) decays mediated by the 6 — > sss transition, one can 
effectively constrain the parameter space of NP, characterised by the strength of the NP coupling and its weak 
phase. 

In this paper we focus upon the relevant decays of type B — > V1V2, which, in this case, are B — > <pK* (with all 
charge combinations) and B s — * </></>. These modes, in particular the former, have been discussed in the context 
of specific NP models |3] as well as in a model-independent way, including possible modifications of low-energy 
QCD dynamics 0J|S]- For B —> (j>K* , data [Hj |7] exists on branching ratios (BR), different CP asymmetries, 
and different polarisation fractions (see Tabic HJ. The error bars are still large but hopefully a much better 
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situation will arise in a few more years. For the B s decay we have only some preliminary data [7] on its BR, 
but LHC-B should do a more thorough job. 

The reason for such an analysis is twofold. First, one can construct more observables than B — » P\Pi or B — > PV 
cases, simply because the final state mesons can be in s, p or cZ-wave combinations. Since the wavefunctions 
have different parity, one can construct CP violating observables even if the strong phase difference between 
various amplitudes be zero. Second, there are a few SM conditions 8 whose violations are relatively simple to 
observe and which will indicate beyond any doubt the presence of NP. 

Experimentally, the anomalous trend persists in B — > V1V2 sector too. The fraction of final states in a lon- 
gitudinally polarised combination is about 50%, whereas one expects this to be dominant over the transverse 
polarisation fractions, which are suppressed by the mass of the decaying quark. At the infinite mass limit, all 
decays should be longitudinally polarised. 

We must emphasise here that the longitudinal polarisation anomaly may turn out entirely to be of SM origin. 
There are discussions in the literature where contributions, neglected so far, have been properly incorporated 
and their effects have been analysed. Needless to say, most of the insights have come a posteriori, after the 
experimental data is announced, but that is only to be expected while dealing with something like low-energy 
QCD. Kagan j^j has shown that the suppression of transverse polarisation is still there in QCD Factorisation, 
but a new strong penguin contribution can lower the longitudinal polarisation fraction in B — > <j)K* . This, 
however, does not affect B — > pK* , where there is no such polarisation anomaly. On the other hand, Beneke 
et al. have shown that there is a significant EW penguin contribution in B — > pK* , and possibly the same 
mechanism works for B — > cf>K* too. Cheng et al. have shown that in perturbative QCD, the longitudinal 
polarisation can go down to 75% if one takes annihilation and nonfactorisable diagrams properly into account 
(without them it is about 92%). All in all, the explanation may lie within the SM, but there is ample motivation 
to look for new physics. 

Here we perform a model-independent analysis of the channel B — > <pK* and extend the analysis to the SU(3)- 
related channel B s — > tf><f>. (Such an analysis, with different set of operators, was also performed in 0, and 
our conclusions are in agreement.) This is the first analysis of all anomalous b — > sss mediated decays in a 
model-independent way, alongwith predictions for B s — > (jxf>. Analyses within the framework of definite models 
(in particular different versions of SUSY) are available, so is a partial model-independent analysis. The problem 
with such partial analyses is that they give different allowed parameter space for the NP for different anomalies. 
Of course, this does not mean that we are evaluating the relevant amplitudes in a model-independent way! 

In fact, this is the second part of the analysis. For the first part involving B — > P\Pi and B — > PV modes, 
we refer the reader to 0. We stress that we had to redo the analysis again since the data changed in the last 
few months, in particular the sin(2/3) anomaly. The results have been summarised in Section IV, but we do not 
repeat the formalism which obviously remains unchanged. 

The data on BR, CP asymmetries, and polarisation fractions are taken as input. For the theoretical input, 
the major uncertainty occurs in the calculation of long-distance contributions. We circumvent the problem 
by a rather conservative approach. The NP effective Hamiltonian is characterised by a real positive coupling 
h, a NP weak phase £ (between and 2tt), and a Lorentz structure for the b — ► sss current-current product. 
All short-distance corrections coming from the running to the NP scale to nib are dumped in h, but just for 
simplicity, we assume the NP operator not to mix with the SM ones. The analysis not only gives the allowed 
region in the h-£ plane, but also predicts the range of different observables. Similar predictions are obtained for 
the B s decay channel. 
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II SM and NP amplitudes 



The amplitude for B(p) — > cj)(kx,£i) + K*(k 2 ,e 2 ) can be written as |T^HT3"| 



M = ae\.e* 2 + — (p.ej)(p.e5) + i-A-e^ap^efef (1) 



mTj mi 



with g = fci — k%, and the CP-conjugate amplitude (for _B) has the obvious form 

M = ae\.e\ + -\-{p.el){jp.e* 2 ) - i-^^™// 6 !"^. ( 2 ) 

771 TTlj^ 

where a, b, and c are in general complex quantities, involving, apart from the short-distance effects, the weak 
and the strong phases. 

In the linear polarisation basis, we write 

1 A (~i \ *r _*r i , \ 1 a i— a — \ (St ^ St 



M(M) = M4>>1*2 ~ ^M A \\( A \\H - £ 2 - (+)^A4A ± ) [ef x e^ 1 j .p (3) 
where /5 is the unit vector along K* in the rest frame of 4>, and 

ef=i*.p, sf = e*-e* L p. (4) 

The amplitudes of eq. © are related with those of eq. JIJ by 

A,, = V2a, A = -ax - ^^b(x 2 - 1), A x = 2V2^^cy / x 2 - 1, (5) 
" m B m B 

(and similarly for the barred variables) where 

k x .k 2 m 2 B - m 2 - mj { . 



mim 2 Im^m 



(6) 



A" 



In the SM, for the decay B — > a and c are real and negative, and b is real and positive (assuming 

negligible strong phases), so that A > (the &-term is mass-suppressed), while Aii,A_l < ^3]. Thus one 
expects 4>± = arg(A±/Ao) w 7r, cf>u = arg(A||/Ao) ~ 7r. However, these expectations may change due to the 
presence of final state interactions (FSI). 

An alternative formulation is in terms of the so-called hclicity basis, where the amplitudes arc written in terms 
of Hq and H±, and 

An ± A_l 

H =A , H±= 11 . (7) 

The decay width is given by 

T = A (l^ol 2 + \H + \ 2 + \HA 2 ) = (|A | 2 + |A„ | 2 + |AJ 2 ) (8) 

where k is the magnitude of the three- momentum of either V\ or V 2 . For the experimental observables, the 
amplitudes are normalised in such a way that |A | 2 + |A||| 2 + |Aj_| 2 = 1, and similarly for the A^s [H] 1 - 

To evaluate the transition amplitudes, we use the conventional factorisation (CF) model ^] E3> with the 
standardised matrix elements as shown below (with q = k\ — k 2 and p = k± + k 2 ): 

?0\ _ _ ,»^„ai,|3 



tub + my 2 



1 This normalisation does not affect the definition of the experimental observables in eq. 1141 . 
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(V 2 \A„\B°) = l (e* £ -l qii \( mB + mv2 )A 
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-* ( (P + k 2 )^ - "° B 2 " v » g/i ) (e*.q)- 



r 

m V2 \ , * N ^2 

2 (e -9) • 

q* ) m B + mv 2 

(Vi\qa^ v q\0) = {^k lv - , 

(V 2 \qa^(l + l5 )b\B ) = 2iT 1 e fiUa/3 E* 2 ''p a k 2 3 



-T-2 {^*2^ m B - m V 2 ) - (4-P) (P + fc 2) M } 



+T 3 ( £ ;.p)(fc lM - ^ ( P + k 2 )\. (9) 

I " l B " l V 2 J 

The form factors, calculated in the light-cone sum rule (LCSR) approach, are taken from P2J. Their values are 
given in Section III. 

This gives, in the SM, the transition amplitude as 



M(B~ -► 4>K*-) = M(B° -> 4>K* a ) 

G F 

= i—i=Um l f } {-(e l p.eK')(mB+rnK')Ai^ K (m^) 

2A^ K '{ m \) g 2V B ^ K '(ml) 

+(^-P)( e K*-P)-f : r - ifF^eK-P"^! 1 \} 

{m B +m K *) v (m B +m K ,) 

x VtbV t * s {a 3 + a 4 + a 5 - -(a 7 + a 9 + a 10 )} (10) 

where the symbols have their usual meaning |15j . However, one may question the validity of the NF approach 
for this decay, and indeed calculations based on QCD factorisation J7| or perturbative QCD indicate a 
discrepancy in the predicted BR by about a factor of 2 at the most (8.71 x 10 -6 in QCD factorisation vis-a-vis 
14.86±3 ^ x 1(T 6 in perturbative QCD for the decay B° -> <j)K*°). To account for this, we have allowed the 
SM amplitude to vary by 40% for a fixed iV c = 3 (this is equivalent to a 96% variation in the BR). Also note 
that this mode, like B — ► <fiK, is not iV c -stable, and the final result may have some quantitative variation for a 
different N c 2 . Apart from this variation, all other effects that may change the predicted BR have been taken 
into account by varying the amplitude. 

One may ask whether it would have been prudent to take QCD factorisation or perturbative QCD as the model 
for the SM dynamics and estimate NP effects. We would like to point out that they are models to calculate 
nonleptonic decay amplitudes, just as CF, and though there is theoretical justification for taking a particular 
model for a particular type of decay, there is no reason to believe that for the decays in question, one is definitely 
better than the others. So we stick to CF; however, it is necessary to take into account the differences among 
various models, at least in a rough way. This will provide for a larger uncertainty in the SM prediction. If wc 
still need NP to explain the data, we may be hopeful about its presence. The question is how one does this. A 
good indicator is the predictions for the BRs. We may hope that, roughly, the differences in BRs are reflected 
in the amplitudes in different models. One may like to have some more error margin, since it is better to be 
conservative than over-ambitious. The same is true for experimental data; in fact, we will soon show that if the 
data is taken at 68% confidence limit (lcr), no possible NP can explain all the anomalies. 

For the decay B s {p) — > <p(ki, Si)<p(k2, s 2 ), the transition amplitude is 

G F 

M(B s ^y(j>4>) = i-j=f (t> m 4> {~2{ei.t 2 ){mB s + m < f,)Af 3 ^' t> (ml) 

AA^{ml) g 2V B ^{ml) 
+ (e 1 .p)(e 2 .p)- ■ \ - le^apefeZp 01 ^-. ■ -} 



2 In the analysis, we stick to N c = 3 since the form factors are evaluated 1101 for this value. 
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x VtbV t * s {a 3 + a 4 + a 5 - ^(a 7 + a 9 + a w )}. (11) 
This is analogous to eq. I)l()[l. with an extra factor of 2 for identical particles in the final state. 
The amplitudes will have contributions from SM as well as from NP. Let us write 



A,, = R ie l0 \A = R 2 e l ° 2 1 A_ L = R 3 e 



l»3 



A,, = i? 4 e^ , A = i? 5 e W5 , A± = ifee" 8 , (12) 

where R4S and 0^s include all SM and NP effects (couplings, weak and strong phases). The 18 variables proposed 
by 8 can be written as 

A 00 = \ {Rl + Rl) , A n = ±{Rl + Rl), A ±± = ± + i? 6 2 ) , 

L ±0 = R 2 R 3 sin(0 2 - 6 3 ) - R 5 R 6 sin(0 5 - 6 ), 
L_l|| = RxR 3 sin(6»i - 6 3 ) - R 4 R 6 sin(0 4 - 6 ), 
L|| Q = i?ii? 2 008(0! - 6 2 ) + i? 4 i? 5 cos(0 4 - 5 ), 

Soo = ^( R l~ R i) > s |||| = 2 (^1 _ ^4) ^ s n = 2 (-^3 _ ^i) > 
S±o = i? 2 i? 3 sin(0 2 - 9 3 ) + R 5 Re sin(0 5 - 6 ), 
S_l|| = i?ii? 3 sin(6»i - 3 ) + i? 4 i? 6 sin(0 4 - 6 ), 
S|| = i?ii? 2 cos(6»! - 2 ) - R4R5 cos(0 4 - 5 ), 

Poo = R2R5 sin(2/3 + 2 - B ), p n = R X R A sin(2/3 + X - 4 ), P±± = -R 3 Re sin(2/3 + 3 - 8 ), 
p ±0 = i? 3 i?5 cos(2/3 + 3 ~ B ) + R2R6 cos(2/3 + 2 - 6 ), 
p_L|l = i? 3 i?4 cos(2/3 + 6» 3 - 4 ) + i?ii? 6 cos(2/3 + X - 6 ), 

Plio = i?ii? 5 sin(2/3 + 0x - 5 ) + R2R4 sin(2/3 + 2 - 4 ). (13) 

Here f3 = arg(V t * d ) is the SM weak phase coming in B° — B° mixing. We assume no NP contribution in this 
mixing. For the B s system, (3 S is close to zero in the SM. However, NP of the type b — » sss may contribute 
to B s — B s mixing. Even then, the contribution of NP in mixing, which is in effect a contamination to /3 S , 
can hardly be worth considering. The reason is this. Only a lower bound on the SM amplitude exists. The 
NP amplitude with such a weak coupling as obtained from the decay fit to B — » P\P 2 or B — > PV modes can 
never compete with the SM amplitude. We find sin(2/3 s ), the effective phase from the B s — B s box, to be never 
greater than 0.1. (Similarly, b — ► ccs channels are hopeless to look for new physics.) 

On the other hand, BaBar and Belle collaborations express their data in terms of eight independent variables 
over which a fit is performed. Apart from Jl = Aqo and /j_ = Aj_j_, they are 



A 

A CP 


fE - fE Rl 
fE + fE 


— R5 

+ RV 






A 1 - 

A CP 


fl-fl Rl 


- R l 






fE + fE Rl 


+ RV 






h 




+■ arg{A\\/A Q )) -- 


= ~ (0i - 2 + 


4 - 5 ) , 




= \(arg{A ± /A ) 


+ arg{A 1 _/A )) 


= 5 (0 3 - 2 + 


' 06 - O5) 




= 2 {«rg(A\\/Ao)- 


-arg(A l{ /A )) = 


= ~ (0i - 2 - 


04 + Ob) , 


A0j. 


= \(arg{A ± /A ) 


-arg(A ± /A )) 


= (^3 - S 2 - 


06 - O5) 



(14) 

where we have used a convention opposite to that used by BaBar, Belle and HFAG to define the first two and 
last two variables of eq. (|14|) . These are the constraints that will go as inputs in our analysis. Note that the 
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set {— (^||,7r — <fi_i, — A(^|| , — A0j_} is identical as far as the angular analysis is concerned. If we entertain the 
possibility of NP, there is no reason to keep our analysis confined to the set with values nearest to the SM 
expectation. 

As in we discuss three different types of effective four-Fermi interactions coming from new physics: 

1. Scalar : C new = h s e l ^ s (s a (ci + c 2 j 5 )s a ) (sp(c 3 + 0475)^/3) , 

2. Vector : C new = h v e^ v (s a 7* t (ci + c 2 j 5 )s a ) (s/37p( c 3 + 0475)6/3) , 

3. Tensor : C new = h t e 1 ^ (s^cr^^i + c 2 j 5 )s a ) (spcr^(c 3 + 0475)6,3) . (15) 

Here a and /3 are colour indices. The couplings h s , v ,t are effective couplings (generically denoted as /ijvp), of 
dimension [M]~~ , that one obtains by integrating out the new physics fields. They are assumed to be real and 
positive and the weak phase information is dumped in the quantities £ s<v ,t (again, generically denoted as £np), 
which can vary in the range 0-2ir. Note that they are effective couplings at the weak scale, which one may 
obtain by incorporating all RG effects to the high-scale values of them. The couplings C1-C4 can take any values 
between —1 and 1; to keep the discussion simple, we will discuss only six limiting cases: 



W (s + 


P) 


X 


(S + P) [(V + A) x (V 


+ A)] 


Cl 


= 1, 


C2 = 


1, C 3 = 


1, c 4 = 1; 


(ii) (S + 


P) 


X 


(S-P) [(V + A) x (V 


-A)} 


Cl 


= 1, 


C2 = 


1, C 3 = 


1, c 4 = -l; 


(Hi) (S — 


P) 


X 


(S + P) [(V — A) x (V 


+ A)} 


Cl 


= 1, 


C2 = 


-1, c 3 


= 1) c 4 = 1; 




P) 


X 


(S-P) [(V-A)x(V 


-A)} 


Cl 


= 1, 


C2 = 


-1, c 3 


= 1, c 4 = - 








(v) (T + PT) x (T ■+ 


-PT) 


Cl 


= 1, 


C2 = 


1, c 3 = 


1, c 4 = 1; 








(vi) (T - PT) x (T - 


-PT) 


Cl 


= 1, 


C2 = 


-1, c 3 


= 1, c 4 = - 



This choice is preferred since the 1 — (+)7s projects out the weak doublet (singlet) quark field. For the doublet 
fields, to maintain gauge invariance, one must have an SU(2) partner interaction, e.g., s(l — 7s)s must be 
accompanied by c(l — 75)0. No such argument holds for the singlet fields. In the above equation, PT denotes 
a pseudotensor structure, characterised by 07^75. 

The tensor current was not considered in Neither this form nor its Fierz-reordered form can contribute to 
B — * <j)K. Now that with the latest data [7j one does not imperatively need NP for the B — ► cf>K sector, one may 
feel justified to include this structure as well. Note that only ct a " / (1 + (— )7s)<8> 07^(1 + (— )7s) structures are of any 
interest; the other two, after Fierz reordering, do not generate any scalar or pseudoscalar currents and hence can 
affect none of the B -> PiP 2 or B -> PV modes. In fact, the four-quark current (T + (-)PT) X (T - (+)PT) 
vanishes, which can be checked from the identity cr A "'75 = —(i/2)e^ va P<j a p, with g 0123 = — 1 (though the 
factorised matrix element need not vanish). 

We have chosen the interaction in a singlet-singlet form under SU(3) C . The reason is simple: one can always make 
a Fierz transformation to the local operator to get the octet-octet structure. Note that the forms (S + (— )P) x 
(S + (-)P) generate tensor currents under Fierz reordering. Such currents were not important in |2] since at 
least one of the final state mesons was a pseudoscalar. Here it will be important since both the final state mesons 
are spin-1 objects, and as we will see, the tensor currents play a crucial role in bringing down the longitudinal 
polarisation fraction of B — > 4>K* . Since no such tensor current is available for (S + (-)P) X (S — (+)P) type 
operators, or the vector-axial vector operators, there is no lowering of the longitudinal polarisation fraction. 

We have kept the strong phase difference between the SM and the NP amplitudes a free parameter. The short- 
distance strong phase, coming from the imaginary parts of the respective Wilson coefficients, are calculable but 
small. The long-distance strong phase, coming mostly from final-state rescattering, is a priori not calculable, 
but since there are not too many final states of identical quark configuration, the strong phase is expected to 
be not too large. However, there should not be any correlation between the strong phase in B — > <f>K* and the 
strong phases in B — > </>K or B — > rf ^K^*\ the channels discussed in but the strong phase of B s — > <p4> can 
be related to that of B — > <\>K* by SU(3) symmetry. We will assume the breaking of flavour SU(3) to be small 
and take equal strong phases in both these cases. The results are not at all sensitive to a precise equality. 
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The a, b, and c terms of the NP amplitudes (eq. JTJ) for the decay processes B — > and i? s — * <^></> take the 
following form: 

£> — > (the contributions are same for neutral and charged channels): 
(i) Scalar-pseudoscalar channel 

a NP = i [-f^m^CtCi - c 2 c 3 )(m B + m K ,)A 1 + 2fJc 1 c i {m 2 B - m 2 K ,)T 2 ] -^-^ NP , 



>NP 



CNP 



f^m^cxCi - c 2 c 3 ) 



2A 2 



t,m^{C2CA - C1C3) 



m B + m K * 
V 



m B + rriK* 



Cl c 4 fJ AT 2 + 4T 3 2 2 



4N r 



(17) 



Note the 1/N C suppression; this channel can only contribute to the decay after reordering. Also, note that for 
(S + (-)P) x (5+ (— )P) structures only the terms with tensor form factors survive, which is obvious since they 
do not generate any vector or axialvector currents. 

(ii) Vector-axial vector channel 



a-NP = Um^ [(cic 4 -I- c 2 c 3 ) + 4iV c cic 4 ] (m B + m K * )A\ 



bNP 



cnp 
m 2 



ll>NP 

4JV ' 



-f^m^ [(cic 4 + c 2 c 3 ) + AN c cic 4 } 
-f^m^ [{c 2 c 4 + C1C3) + 4A^ c ciC 3 ] 



2A 2 



(tub + rriK* ) AN C 



y_ h NP r i £NP 



(tub + m K , ) AN C 



(18) 



(hi) Tensor-pseudotensor channel 



aNP 
bNP 



m 



B 



CNP 



-i [ 2 f<p c l c i( m B - m K*) T 2] ( 1 



Cl c 4 /J 4T 2 +4T; 



m- 



3 — 2~~ 

TOT, — TO 



2N r 



2N r 



h N pe iiNP , 



K * 



hN P e iiNp . 



2N r 



tiNpe 



(19) 



The expressions for B s — > (jxfi are analogous, with the obvious replacements B — > _B S , ii"* 
factor of 2; see, for comparison, eqs. (fTT)f> and fTTjl. We do not tabulate them separately. 



and an extra 



III Theoretical and Experimental Inputs 

The experimental data, taken from |7j, is shown in Tabled The numbers are quoted for B — > <f>K* (neutral 
mode) while the corresponding numbers for charged B decay, wherever they exist, are given in parenthesis. The 
error margins are shown at la confidence limit (CL), while for the analysis, we have taken a more conservative 
approach and kept the error margins at 2a. We do not use the numbers that are derived from the primary 
measurements assuming the validity of the SM, mostly A|| n , L|| , and various £s. Note that since we are interested 
only in the b — ► sss transition, no other decay modes (like B — > pp) have been taken into consideration. 



While the BR(£? — > 4>K*) is in the expected ballpark, BR(_B S — > (jxf>) is smaller than expected. The amplitude 
for the latter is twice that of the former (identical particles in the final state) times the SU(3) breaking effects, 
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which gives an enhancement of the BR of the latter by roughly a factor of 5.5-6. But the number of events for 
B s — > (jxf) is small; it is compatible with zero at 95% CL! For analysis, we take this particular data with a 3tr 
error bar on the higher side instead of the usual 2a, just to be cautious over the preliminary nature of the data. 



Observable 


Value 


Observable 


Value 


Br(B -> cj)K*) 
Aoo 

h 

A0|| 


(9.5 ±0.9) x 10" 6 
((9.7 ±1.5) x 10~ 6 ) 
0.48 ±0.04 
(0.50 ±0.07) 
2 3fi+ 018 

z -°°-0.16 

(2.10 ±0.28) 
-0.01 ±0.08 
-0.03 ±0.18 


Bv(B s -> #) 
f± = A±± 

A 1 - 

A CP 


(141?)) x 10" e 

0.26 ±0.04 
(0.19 ±0.08) 

2.49 ±0.18 
(2.31 ±0.31) 

0.16 ±0.15 
-0.03 ±0.18 



Table 1: Data on B — > <fiK* modes, from |7j. Our convention of defining CP asymmetries is opposite to that of 
HFAG, see text. We do not show, for obvious reasons, those observables which are not directly measured but 
estimated using the validity of the SM. 

Apart from the BRs and CP asymmetries, we also use the following results from [7J: 

• sin(2/3) from charmonium modes: 0.685 ± 0.032; 

• sin(2/3) from B — > Ks<f> transitions: 0.47 ± 0.19 (the results do not show any qualitative change if we use 
the combined b — > sss result: 0.50 ± 0.06; however, the averaging is a bit naive [7J and should be used 
with caution); 

• 43.8° < 7 < 73.5° at 95% CL ^2]; this is needed for a reevaluation of the constraints on the allowed 
parameter space (APS) of new physics as found in 

The CKM elements Vt s and Vtb are taken from |19| . with only the unitarity constraint imposed. 

The constituent quark masses, in GeV, are taken to be ^H] 

m u — ma = 0.2, m s = 0.5, m c = 1.5, m& = 4.88, (20) 

though the final result is totally insensitive to the precise values. The constituent quark masses are independent 
of the renormalisation scale. The Wilson coefficients, evaluated at the regularisation scale \i = m&/2 w 2.5 
GeV, are also taken from |15| . The corresponding current quark masses, which appear in the Dirac equation 
for quarks while evaluating the hadronic matrix elements, are (in GeV) 

rn u = 0.0042, m d = 0.0076, m s = 0.122, m c = 1.5, m b = 4.88 (21) 

at n = 2.5 GeV (see table I of the first paper of JS]). The light quark masses will obviously shift if we evaluate 
them at some other scale, say 1 GeV, and they also depend on how we evaluate them (e.g., lattice QCD or 
QCD sum rules). The decay constant of <f>, through vector and tensor currents, are defined as 

(0|g^»|^(p,A)) = ifj (eWp, - eS A V) , (22) 

and their numerical values (in GeV) are ^B] 

U = 0.231 ± 0.004, /J = 0.200 ± 0.010. (23) 
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The form factors, evaluated in the light-cone sum rule (LCSR) approach, are JHI 

B->K* : Ai = 0.292, A 2 = 0.259, V = 0.411, T X =T 2 = 0.333, T a = 0.202, 
B s ^(j): Ax = 0.313, A 2 = 0.234, V = 0.434, T x = T 2 = 0.349, T 3 = 0.175. (24) 

These form factors include the full twist- 2 and twist-3 and the leading order twist-4 contributions. These 
numbers are for q 2 — 0. For nonzero q 2 , they change by about 10%. We take the q 2 — values for our 
numerical evaluation. They are evaluated for iV c = 3. 

IV Results 



Lorentz 
structure 


(xlO" 8 ) 


t,NP 
(°) 


{S-P)x {S + P) 


0.0-3.0 
(0.5 - 3.0) 


180 - 360 
(185 - 340) 


(S + P)x (S + P) 


0.0-4.5 


0-180 


{S-P)x (S-P) 


0.0-4.5 


180 - 360 


{S + P)x (S-P) 


0.0-3.0 


0-180 


(V - A) x (V + A) 


0.0-0.8 


180 - 360 


(V + A) x (V + A) 


0.0-2.1 
(0.25-0.7) 


0-180 
(30 - 140) 


(V — A) x (V — A) 


0.0-0.6 


180 - 360 


(V + A) x (V - A) 


0.0-2.2 
(0.5 - 0.8) 


0-180 
(40 - 130) 



Table 2: Allowed parameter space from the analysis of B — ► P\P 2 and B — ► PV modes, mediated by b — ► sss 
transition (upgrade of |2] in view of the Summer 2005 data [7]). All error bars are taken to be at 2a. With la 
error bars, only three structures survive, as shown in parenthesis. 

Before we embark on an analysis of B — > V\V 2 modes, let us revisit the results of 2 in the light of Summer 
2005 data. The main change is the prediction of sin(2/3) from B — > c/)Ks' this is now less than 2a away from 
the charmonium result. Thus, a nonzero NP amplitude based on this data alone is no longer necessary (and 
hence tensor currents are a possibility). Of course, the branching ratios of the ?/ modes are still too large, 

but there is no way one can reconcile that with a pure NP amplitude; there must be some dynamics beyond the 
naive valence quark model (Lipkin has suggested that this is due to interferences between B — * Kr/s 
and B —> Kr\\ amplitudes, constructive for rj and destructive for r\. Recently, an analysis based on the Soft 
Collinear Effective Theory (SCET) claims this to be supported in that framework However, SCET is yet 
to provide numerical data on all the modes taken in our analysis.) The allowed regions for different Lorentz 
structures is shown in Table 2. Note that h^p can be vanishingly small. However, the upper limits, which are 
controlled by the BRs, remain unaltered. There does not seem to be any pressing need to introduce new physics 
from this data alone, but we will soon find that the longitudinal polarisation anomaly forces us to consider the 
NP option seriously. 

In Table [3 we show the principal allowed parameter space (APS) for different Lorentz structures. For each 
structure, there is a subdominant parameter space with very small h^p, but £,np in the opposite half-plane 
(i.e., £np + 7T modulo 2ir). For example, for the structure (S + P) x (S + P) there is an APS with very small 
h>NP ~ 10~ 9 and 7r < £np < 2tt. This was absent in [2], but now that sin(2/3) from B — > 4>Ks at 2a may 
overshoot the charmonium value, there is a scope for opposite interference. This is a general trend for all 
structures. However, these regions do not survive the B — > ViV 2 analysis. Also note that with la error bars, 
three structures survive, in contrast to only one as found in [2]. Again, they disappear after the B — > ViV 2 
analysis. 
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Set 
No. 


Lorentz 

□ 11 1.1 V.. 1 : l.l L \j 


h-NP 
l Alu ) 


\ ) 


\ ) 


I 


(S + P)x(S + P) 


2.0 - 3.6 


- 30 
45 - 65 
150 - 180 


90 - 145 
25 - 35 
270 - 325 


II 


(S-P)x(S- P) 


2.8 - 3.5 


271 - 278 

ooo ouu 


216 - 223 


III 


(T + PT) x (T + PT) 


0.14-0.36 


- 25 

- do- 

1 50 - 1 80 

- do- 
45 - 70 


140 - 170 
270 - 320 
90 - 140 
320 - 350 
205 - 215 


IV 


(T - PT) x (T — PT) 


0.20 - 0.27 


- 15 
90 - 100 
158 - 180 


125 - 135 
215 - 225 
305 - 320 



Table 3: Allowed parameter space for different Lorentz structures of new physics. 



With the expressions for BRs and CP asymmetries at hand, we perform a scan on the new physics (NP) 
parameters hnp and £,np- The starting ranges for each Lorentz structure are shown in Table [21 We perform a 
full scan over the tensor structures. This time, however, we must introduce a nonzero strong phase difference 
between the SM and the NP amplitudes (for simplicity, we assume this difference to be the same for all angular 
momentum channels). The reason is that both <f>n and (f>j_ differ from the SM expectation of ir, even taking the 
uncertainties A(f> into account. It is easy to see that if the strong phase difference Spfp is zero (modulo 7r), both 
(pu and 4>±_ retain their SM expectation values. If the strong phase is generated from rescattering, it should 
not be related in any way from the strong phases in B — > P1P2 or B — > PV channels; but from SU(3) flavour 
symmetry, we expect the same Snp (at least to the leading order) for both B — > <f>K* and B s — > (f>cj). In our 
analysis, we take them to be the same. This effectively reduces our parameter set to /ijvp, £/vp and Snp- 

Using the inputs discussed earlier, we find the APS for these three parameters for different Lorentz structures. 
The result is shown in Table 01 Note that only two scalar-pseudoscalar and two tensor-pseudotensor channels 
survive. This is due to the fact that only these channels, under Fierz reordering, generate a tensor current, 
which helps to bring down the longitudinal polarisation fraction Aoo- This is in agreement with [5]. These 
results are obtained with 2a error bars; nothing survives at la CL. 

In Tables 01 and we show our main results: our expectations for the four allowed sets. The results are quite 
similar for B — > <f>K* and B s — > <jxf) 1 which is nothing but the manifestation of a rough SU(3) symmetry. In 
particular, we expect a similar suppression of Aoo for B s — > (jxj) too. 

At this point, one may like to look at the SM predictions for these observables. In the conventional factorisation 
scheme, they are: 

(i) B -> (j>K*\ 

Aoo = 0.893, Aj_j_ = 0.051, K n = 0.056, Ap = 0.446, 
poo = -0.035, p ±± = -0.612, = -0.038, 

p±o = 0.312, pyo = 0.306, p X || = 0.078. (25) 

(ii) B s -> 00: 

A 00 = 0.948, A x _l = 0.009, A u = 0.044, A|| = 0.406, 

p±o = 0.184, pj.11 = 0.039. (26) 
All the other variables are zero. Again, this is in the conventional factorisation model; nonfactorisable and 
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Observable 


Set I 


Set II 


Set III 


Set IV 


A 00 

A±± 

A llll 
A± 
Ayo 
A j. || 


0.40 -> 0.56 
0.21 -> 0.29 
0.23 -» 0.32 
-0.15 -> 0.25 
0.25 -» 0.71 
0.0 


0.40 -» 0.56 
0.18 -» 0.26 
0.21 -> 0.40 

-0.06^ 0.16 
0.46 -» 0.70 

-0.30 -> 0.10 


0.40 -» 0.56 
0.21 -» 0.29 
0.23^ 0.32 
-0.25 -> 0.25 
0.25 -> 0.70 
0.0 


0.41 -> 0.56 
0.18^ 0.27 
0.23^ 0.41 

-0.06 -» 0.18 
0.46^ 0.75 

-0.29 -> 0.18 


^00 

s±± 

s llll 

£±o 

S|| 

Si II 
- 1 - ll 


-0.07^ 0.03 
-0.01 -» 0.03 
-0.01 -> 0.04 
0.22 -» 0.64 
-0.21 -> 0.22 
0.0 


-0.04 -» -0.005 
0.01 -» 0.14 

-0.09 -> -0.006 
0.44 -» 0.60 

-0.37^-0.05 
0.005 -» 0.25 


-0.07^ 0.03 
-0.01 0.03 
-0.01 -» 0.04 
0.24^ 0.62 
-0.20^ 0.22 
0.0 


-0.04^ 0.12 
-0.03^ 0.13 
-0.09^ 0.18 
0.37^ 0.60 
-0.38^ 0.12 
0.005 -> 0.23 


Poo 

P-L-L 

Pllll 
P±0 

P||0 

p±\\ 


-0.25 -> 0.05 
-0.41 -» -0.18 
-0.27^ 0.05 

0.42 -» 0.70 
-0.10^ 0.42 

0.18 -> 0.70 


-0.25 -» 0.05 
-0.50 -» -0.25 
-0.21 -» -0.04 
0.47^ 0.67 
0.16 -» 0.27 
0.34 -> 0.60 


-0.3 -> 0.05 
-0.40^ -0.18 
-0.35 -» 0.07 

0.40 -» 0.70 
-0.10^ 0.50 

0.20^ 0.65 


-0.26 -» 0.04 
-0.50^ -0.25 
-0.35 -» 0.0 
0.46^ 0.67 
0.14 0.42 
0.34^ 0.58 



Table 4: Observables for B -» 0K* . 



Obscrvablcs 


Set I 


Set II 


Set III 


Set IV 


A o 
A±± 

A llll 
A± 

A l|0 

A±|| 


0.39 -> 0.66 
0.14^ 0.32 
0.20 -> 0.33 

-0.70 -» 0.60 
0.10 -> 0.72 

-0.08 -» 0.07 


0.40 -» 0.60 
0.19 -» 0.29 
0.21 -» 0.36 

-0.06^ 0.16 
0.52 -> 0.70 

-0.20^ 0.04 


0.20^ 0.65 
0.15 -» 0.42 
0.20^ 0.38 

-0.70 -» 0.60 
0.10^ 0.73 

-0.08 -» 0.07 


0.38 -» 0.58 
0.18 -» 0.31 
0.22 -> 0.36 

-0.09 -» 0.18 
0.49 -> 0.74 

-0.20 -» 0.06 


Soo 

£±± 

s llll 
£±o 

s ±ll 


-0.026 ^ 0.015 
-0.025 ^ 0.01 
-0.015^ 0.04 
-0.24^0.70 
-0.25 -> 0.25 
-0.06 -» 0.08 


-0.04 -» -0.002 
0.01 -» 0.09 

-0.06^ -0.03 
0.44 -» 0.66 

-0.37^ -0.04 
0.004^ 0.16 


-0.030^ 0.013 
-0.024 -> 0.01 
-0.013^ 0.04 
-0.22 -> 0.70 
-0.30 -» 0.25 
-0.07^0.08 


-0.04^ 0.008 
-0.018^ 0.09 
-0.06 -» 0.01 
0.36 -» 0.65 
-0.38 0.12 
0.004^ 0.15 


Poo 

P_L_L 
p\\\\ 

P±o 

P||0 

PM\ 


-0.26 -» 0.25 
-0.05 -> 0.11 
-0.22 -> 0.30 
0.28 -> 0.62 
-0.46 -» 0.25 
-0.30 -> 0.50 


-0.13 -» 0.15 
-0.15 -» 0.06 
0.06 -» 0.22 
0.40 -» 0.75 
-0.40 -» -0.22 
0.30 -» 0.55 


-0.26 -» 0.22 
-0.04 -> 0.11 
-0.20 -> 0.24 
0.26^ 0.74 
-0.47^0.30 
-0.30 -» 0.80 


-0.15 -» 0.18 
-0.15 -> 0.05 
-0.14 -> 0.25 

0.37^ 0.77 
-0.42 -> -0.14 

0.28 -> 0.60 


A CP 

Acp 
^ 
<t>\\ 
5cf>± 
64> n 


-0.15 -» 0.10 
-0.2 0.1 
-1.2 -> 0.8 
-1.1 -> 0.6 
-1.5 -» 1.5 
-1.2 -» 0.5 


-0.22 -» -0.02 

0.03 -» 0.25 
-0.95 -» -0.82 
-0.97^ -0.55 
-0.15 -» 0.35 
0.08 -» 0.36 


-0.20 -» 0.10 
-0.20 -» 0.10 
-1.2 -> 0.8 
-1.2 -» 0.6 
-1.3 -» 1.3 
-1.2 -> 0.4 


-0.23 -» -0.02 

0.03 0.25 
-0.98 -» -0.65 
-0.96 -0.46 
-0.25 -» 0.40 
-0.18 -» 0.3 



Table 5: Observables for B s — > 00. 
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annihilation contributions may reduce A o of B — > cf>K* to about 0.75 More variables are zero for the 

latter case since we have two identical vector mesons and the B s — B s mixing phase is close to zero in the SM. 

From tables^and|31 it appears that possible four structures can be divided into two sets: one with (S+P) x (S+ 
P) and (T+PT) x (T+PT), and the other with (S-P) x (S-P) and (T-PT) x (T-PT). Precise measurement of 
all the observables should be able to discriminate between the two sets, but considering the respective numbers, 
this is a more than formidable job. (An almost impossible task is to discriminate between the pair of a given set. 
One way out is to look for anomalies in semileptonic decays and analyse the angular distribution of the emitted 
leptons.) On the other and, nonzero values of most of these observables will point to new physics. Note that the 
values of 4>\\ and <p± are modulo 2tt, and the ambiguity of {<j>\\, A</>ii, A</>^} <-> {— (fin, tt — cf>j_, — A0m, — A<^} 
is still there. 

This analysis makes some of the more favourite models of new physics less so. A prime example is R-parity 
violating supersymmetry, which generates only (S + (-)P) x (S — (+)P) type interactions, but not those that 
survive our analysis. The NP particles may be directly detected at the LHC if the corresponding dimensionless 
couplings of the full theory are perturbative. For example, if it is ~ 0.1, then we expect new particles to be 
about 200-400 GeV, perfectly in the range of LHC. This is true even for tensor currents if the tensor structure 
appears from an underlying radiative effect with loop suppressions coming into play. 

We would still like to mention again that the analysis should not be taken as as irrefutable proof for and against 
certain types of new physics models. We have tried to constrain the parameter space for generic NP models, 
but unfortunately the SM uncertainty is still inordinately large. We have taken a middle-of-the-way approach 
and tried to take into account the variations in predictions of different models (for amplitude calculation) by 
incorporating some uncertainties in the SM prediction by hand. The uncertainties are so chosen as to cover the 
BR predictions of different models. Hopefully, the estimates on NP parameter space err on the conservative 
side. What one needs to make the arguments watertight is to have more control on the theoretical uncertainties. 
It is imperative that we have, as soon as possible, a particular theoretical model whose results can be relied 
upon upto a certain order. This model should have a sound theoretical justification and should not concentrate 
on reproducing the experimental results alone (in particular those that may contain hints of new physics). With 
this, and more data, can one hope to refine this analysis. 
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